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Adopting the framework of the Jaynes-Cummings model with an external quantum field, we obtain 
exact analytical expressions of the normally ordered moments for any kind of cavity and driving 
fields. Such analytical results are expressed in the integral form, with their integrands having a 
commom term that describes the product of the Glauber-Sudarshan quasiprobabihty distribution 
functions for each field, and a kernel responsible for the entanglement. Considering a specific 
initial state of the tripartite system, the normally ordered moments are then applied to investigate 
not only the squeezing effect and the nonlocal correlation measure based on the total variance 
of a pair of Einstein-Podolsky-Rosen type operators for continuous variable systems, but also the 
Shchukin-Vogel criterion. This kind of numerical investigation constitutes the first quantitative 
characterization of the entanglement properties for the driven Jaynes-Cummings model. 



> 

oo 
o 



keywords: Driven Jaynes-Cummings Model, Entanglement, Inseparability Criteria 
PACS: 03.65.-w; 03.65.Ud; 03.67.Mn 



I. INTRODUCTION 



o 

\Q • Entanglement and nonlocal correlations are abstract concepts that naturally appear in quantum mechanics when 

the superposition principle is applied to composite systems. Nowadays, these concepts play an essential role not 

"p^ . only in the quantum computation scenario [If and quantum information theory .2|, but also in the context of rel- 
O |' ativity theory [3]. For instance, one of the main tasks of quantum information theory (QIT) is to develop a quan- 
titative characterization of the entanglement properties and quantum correlation effects for multipartite physical 
systems described by continuous and/or discrete variables [1,0]. In this sense, recent theoretical and experimental 
developments in QIT, which are based on the continuous variable regime and with emphasis on quantum opti- 
cal implementations involving quadrature amp litudes of the electromagnetic field, have appeared in the literature 
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. From the theoretical point of view, Duan et ai. [6| proposed a simple inseparability criterion based on the calculation 
of the total variance of a pair of Einstein-Podolsky-Rosen (EPR) type operators for continuous variable states. Basi- 
cally, this criterion states that 'for any separable continuous variable states, the total variance is bounded from below 



by a certain value resulting from the uncertainty relation, whereas for entangled states this bound can be exceeded.' 
Consequently, the violation of this bound can be interpreted as 'a sufficient condition for inseparability of the states.' 
Similarly, Simon Q adopted the Peres-Horodecki criterion of positivity under partial transpose in the context of 
separability of bipartite continuous variable states, and showed that this mathematical operation admits a geometric 
interpretation as mirror reflection in phase space. Thus, exploiting the consequences of such geometric interpretation, 
the author derived 'uncertainty principles stronger than the traditional ones, to be obeyed by all separable states.' 
Following the same mathematical approach adopted in Giedke et al. [l2| showed that two-mode squeezed states 
maximize EPR-like correlations for a fixed amount of entanglement; besides, this result was used to determine the 
entanglement of formation for all symmetric gaussian states describing bipartite systems. More recently, Shchukin 
and Vogel [H| established an important set of inseparability criteria for bipartite quantum states which is formulated 
in terms of observable moments associated with a variety of quantum states. To this end, the authors derived a 
hierarchy of necessary and sufficient conditions for the negativity of the partial transposition of bipartite quantum 
states (or sufficient conditions for entanglement) that generalizes some previous purposes from literature. 
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Now, from the experimental point of view, Josse et al. [14| produced in laboratory both quadrature and polarization 
entanglement via the interaction between a linearly polarized coherent field and a cloud of cold cesium atoms placed 
in a high finesse optical cavity. To demonstrate continuous variable entanglement in this beautiful experiment, the 
authors used the theoretical tools established in [(| [l2[ - namely, the inseparability criterion and the entanglement of 
formation - with great success in the experimental data analysis. Furthermore, this experiment opened new windows 
of investigation in similar physical systems which describe the matter- field interactions. A feasible physical system 
to generate continuous variable entanglement is given by the Jaynes-Cummings model (JCM), such theoretical tool 
being typically realized in cavity-QED experiments involving Rydberg atoms crossing superconducting cavities in 
different frequency regimes and configurations, with relaxation rates small and well understood [23[. Beyond these 
fundamental features, it is worth mentioning that multipartite gaussian states have potential applications for quantum 
teleportation [H, [25[ and quantum criptography [2(| . 

Many authors have investigated the two-mode and driven JCM in different theoretical contexts and predicted new 
interesting results (for example, see [13, [H, [2t| [3(J HH [H, [H, [H, [H, HE H3) HI, HI) 53] )■ However, a quantitative 
characterization of the entanglement properties and a detailed analysis of the nonlocal correlation effects among 
the constituents (fields and atoms) of these similar physical systems have not frequently appeared in the literature, 
and represent two important additional tools in the comprehension of the atom-photon interactions which deserve 
to be carefully studied. Thus, the main aim of this work is to present the first quantitative characterization of the 
entanglement properties for the driven JCM, which is based upon the inseparability criteria for continuous variable 
systems developed in [y, [16[. To this end, we assume the atom is initially prepared in the excited state, and the cavity 
and external fields are in the diagonal representation of coherent states. Next, we use the mathematical procedure 
developed in [4fj| to obtain exact analytical expressions of the normally ordered moments for any kind of cavity and 
driving fields. These analytical results are then expressed at the integral form with their integrands presenting a 
common term that describes the product of the Glauber- Sudarshan quasiprobability distribution functions for each 
field, and a kernel responsible for the nonlocal correlations among the constituents of the tripartite system. To 
illustrate our results we fix both the cavity and driving fields in a coherent state. In particular, we attain new 
results within which some of them deserve to be mentioned: (i) we establish a link between squeezing effect and 
entanglement for different values of detuning frequency in the driven JCM; (ii) we show that EPR uncertainty can not 
be considered the first quantitative characterization of the entanglement properties of the system under investigation; 
and finally, (iii) we present a numerical evidence that supports the hierarchy of necessary and sufficient conditions for 
inseparability derived by Shchukin and Vogel [16j |. 

This work is organized as follows. Section II describes the mathematical procedure used in the derivation of exact 
analytical expressions of the normally ordered moments for the driven JCM. In Section III we apply our results in the 
analysis of the squeezing effect for each field fixing, for convenience, both the cavity and driving fields in a coherent 
state. Moreover, we also analyse (by means of a numerical investigation) some recent proposals of inseparability 
criteria for physical systems described by continuous variables (in particular, we consider the EPR uncertainty and 
the Shchukin- Vogel criterion). Section IV contains our conclusions. Finally, Appendix A shows the main steps to 
calculate the exact analytical expressions for the generalized moments. 



II. DERIVATION OF THE NORMALLY ORDERED MOMENTS 

In many textbooks on quantum optics [4l| , the normally ordered moments are generally defined in terms of an 
auxiliary function (also denominated as the normal characteristic function) which describes the normal ordering of 
creation and annihilation operators of the electromagnetic field - i.e., 



A^(£,C;^Tr[p(t)e« ct e-«* c 



p(t) being the density operator of the system under investigation. The connection between normally ordered moments 
and normal characteristic function is promptly established through the relation 
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for {r, s} G N. However, depending on the circumstances associated with a particular physical system, the connection 
between normally ordered moments and Wigner characteristic function is more appropriate to our needs, 



<ct^)c s (i)> = (-1)' 



ir+s 



9C9£* 



■e*l«l'AW(€,r;t) 



(2) 



FIG. 1: Experimental scheme used in the description of the driven Jaynes-Cummings model. It consists of a high-Q supercon- 
ducting cavity containing two electromagnetic fields prepared in the coherent states, and a beam of excited atoms interacting 
nonresonantly with both the cavity field and the external quantum (driving) field. 



where A^(£, = Tr [p(i)D c (£)] describes the symmetric ordering of the creation and annihilation operators, and 
D c (£) = exp(£ct — £*c) is the displacement operator. Thus, the initial task for both situations consists basically in the 
calculation of characteristic functions associated with a physical system of interest described by the density operator 
p(t). Consequently, the exact analytical expression of the normally ordered moments will depend on the derivatives 

of A { N c) (Z,C;t) or e£KI 3 A$ (£,£*;*) with respect to £ and <f at the point £,<f = 0. 

Now, let us consider a feasible physical system modelled by a two-level atom interacting nonresonantly with a 
single-mode cavity field, and driven additionally by an external quantum field through one open side of the cavity 
(the proposed experimental apparatus is sketched in Figure 1). Within the dipole and rotating- wave ap prox imations, 
the dynamics of the atom-cavity system is governed basically by the Hamiltonian H = Ho + V, where [40( 

Ho = huS + -hajer z , 

V = ~M<r z + hK ett (A^(T- + A<r + ) . 

Here, S = n a + rib is a conserved quantity (n a = a^a and rib = b b correspond to the photon-number operators of 
the cavity and driving fields), A = e a a + £bb represents a quasi-mode operator with e a (b) = K a (b)/ K eit ( K c ff = K a + K b 
is an effective coupling constant between the atom and an effective field described by the quasi-mode operator), and 
5 = ujq — uj is the detuning frequency between the atomic transition frequency luq and the cavity field frequency u> (in 
particular, we assume the resonance condition between the cavity and driving fields). In addition, the atomic spin- flip 
operators <r± and er z are defined as cr + = \e){g\, <j = \g)(e\ and <r z = |e)(e| — |<?)(g| (|<?) and |e) correspond to the 
ground and excited states of the atom), which obey the commutation relations [<r z , cr±] = ±2cr± and [<x+, er_] = er z . 
Once that [Ho,V] = 0, we obtain, in the interaction picture, the Hamiltonian Hj nt = V, and this fact allows us to 
describe the system through the well-known nonresonant JCM Hamiltonian for an atom interacting with the quasi- 
mode A, and whose coupling constant is given by K off . Consequently, the unitary time-evolution operator is the usual 
nonresonant JCM time-evolution operator, namely U(i) = exp(— iVt/h). 

After this brief introduction of the driven JCM, let us describe the density operator p(t) = U(i)p(0)U^(i), p(0) 
being the density operator of the system at time t — 0. For this purpose, we assume the atom is initially prepared in 
the excited state and the cavity (a) and driving (b) fields are in the diagonal representation of coherent states - i.e., 
the initial density operator is written as p(0) = p at (0) <8> p a b(0), with p at (0) = |e)(e| and 

Pab(0) = [[ — -^o — ~ -Pa(a a )-Pb(ab)|cka, Ob)(a a , Oh | , 



where P(a) represents the Glauber-Sudarshan quasiprobability distribution function for each field. Thus, the matrix 
elements (t) of the tripartite system in the atomic basis can be evaluated through the mathematical expressions 

p n (t) = U n (i)P a b(0)uU) , 

Pl2 (t) = U n (i)p ab (0)l4(t) , 

p 21 (t) = \J 21 (t) Pah (0)V\ 1 (t) , 

p 22 (t) - U 21 (i) Pab (0)l4(t) . 
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In these equalities, the matrix elements Uy (t) are given by [4l[ 

, . / / — \ 5 sin (t\/3. ) 

Un(t) = oos(tVft")-i- ^ , 

U u (t) 

U 2 i(t) 



1 ^eff A 



U 22 (t) = cos(t^7)+i- 




with <p a 



,.N a 



(<5/2) 2 l, (3 A = (p A + K 2 ff l, and N A = A^A. Note that p(t) describes the exact solution of the 



Schrodinger equation in the interaction picture with the nonresonant driven-JCM Hamiltonian. Using this solution 
we can establish exact analytical expressions for the time evolution of some quantities (e.g., the Wigner characteristic 
function and the normally ordered moments for each field) which allow us to shed more light on the quantitative 
aspects of entanglement in the system under consideration. 

The Wigner characteristic functions for both the cavity and driving fields can be formally expressed in the integral 
form, with their integrands having a common term that describes the product of the Glauber-Sudarshan quasiprob- 
ability distribution functions for each field, and a kernel responsible for the entanglement among the constituents of 
the system - namely 



where 



and 



A#&r;t)= // 2 Qb P a K)P b K)Kg,(a a ,a b ;t) 



(a a , a h ; t) = f/D^, (a a , a b ; t) + 21 (a a , a b ; t) 



(a a ,ab|UJ(t)D c (f)U y (t)| 



(3) 



(4) 



(5) 



The exact analytical expressions of the normally ordered moments for any kind of cavity and driving fields are 
then properly obtained by means of Eq. ([2|) when c = a,b. Henceforth the initial states of the cavity and driving 
fields will be fixed in the coherent states {|^ a (b))} throughout this paper. Such mathematical procedure diminishes 
considerably the technical difficulties found in the double integrals over the complex a a - and ab-planes, because (i) 
the Glauber-Sudarshan quasiprobability distribution functions for each field coincide with P a (a&) = Tr8^(a a — u a ) 
and Pb(ab) = Trd^ 2 \ab — (5^- 2 \z) denotes the two-dimensional delta function) and consequently, (ii) the Wigner 



characteristic function AV (£, £ 



{z-^lc c*-f~j i s equivalent to Kl a ;?(V a , t). Therefore, the normally ordered moments are 



expressed within this context as follows: 



<at'(tX(i)> = (-1)' 



<bf(i)b*(i)> = (-1)* 
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The next task will consist basically in the calculations of the mean values (O for the cavity and driving fields, and in 
the derivation of theirs respective normally ordered moments. 



A. Cavity field 



The mathematical procedure adopted in the derivation process of the mean values ([5]) for the cavity field is based 
upon the well-established results for su(2) Lie algebra and its respective decomposition formulas [42j. After lengthy 
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calculations, the analytical expressions for the mean values assume the exact forms [4C 

oo 

m.m'— 

oo 

r (m+l,m' + l) , 
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21 



where the complex functions yCj^ TO /(^ a , Vh) and Y^™^ m '(z'a, t'b) are given by 



ij C m,m'(^a, Vb) = —, 6Xp ( - |e a I/ a + e b ^b| J ka^a + ^b^b 



\2[j—i-\-m) 



and 



Ygr'Va, H.) = exp j -i|C| 2 + 2ilm [e b (e b „ a - e a v h )* £] \ [e a (e a v a + e h v h )* <] m '~ m ' ^ (e^l 2 ) 



In these expressions, L„ (z) denotes the associated Laguerre polynomials, 



F m (t) = cos 



A 7n t 



— l — — sin 

Am 



fi, 



and Gm(i) = — i — — sin 



A m t 



correspond to functions responsible for the time evolution of the mean values, with £1 m — 2 K^ff^/ 771 -j- 1 being the 
effective Rabi frequency and = 8 2 + O^. 
Now, let us introduce the auxiliar function 



which has the following exact expression 



(8) 
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with 



^m'm'^a, ^b) = (^b^a ~ Ca^b)^ (^b^a - Ea^b)]' 
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Here, p = Min[777, r, s — (m' — m)] yields the smallest positive integer element of the set {m, r, s — {m! — rre)}. One 
virtue of this auxiliar function is that it allows us to derive an exact analytical expression of the normally ordered 
moments for the cavity field, i.e., 

I 

(at^)a^)), = J2 ^ t) 

rn—0 

£—1 £ 
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where 

+ v /(m + l)(m' + l) 2iC + i,w+iK, H,) 6^i >ro+1 (^ a , i*) G m (i)G r * n , ( t ), 

.*m&K.^;*) = uCK,n,)6^(".^)^(t)C(<) 

+ v /(m + l)(m' + l) aiC^+l.m+lK, "»0 ©m+l.m'+lK, ^b) G ro (i)G^,(t). 

Note that the infinite sums present in y D^, (i/ a , i) were substituted by finite sums, I being the maximum value 
which guarantees the convergence of this expression (we have fixed I = 60 in the numerical investigations). Indeed, 
exact expressions are reached only in the formal limit I — ► oo. Following, we will establish similar results for the 
external driving field. 



B. External driving field 



Adopting an analogous mathematical procedure for the external driving field, we verify that jj £* (y a , i-^; t) 
have similar structures to the previous cases but differ in the dependence on the variables £ and £* - i.e., the complex 

function Y^™»™ 1 (v a , v^) should be adequately replaced by 



rfira.rn') , \ 
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for each i,j = 1,2. Consequently, the auxiliar function 

®^K^) = (-ir 
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has a central role in the present approach since its analytical expression 
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- Efc!4 ro - fc) (0)4 r - fc )(0)i^-l:r fe) (0)3£WK,^), 

2 A' 



fe=0 



e a |e a ^b - £b^ a | 



allows us to derive exact results for (h^ r (t)h s (i)) through the Wigner characteristic function A^ (£, It is worth 

mentioning that the validity of this auxiliar function obeys the same conditions verified for Eq. |[8j . 

Now, after some tedious calculations, the exact analytical expression of the normally ordered moments for the 
driving field can be expressed as follows: 

e 

(tf r (t)b s (t)U = ^3^K^b) b H^K,^;f) 

e-i e 

m'—0 m— m'-f-l 

e-i £ 

+ E E 3^K,^)b*^K,^;t), (ii) 
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where 

bSfc&K, ^ t) = nC£ m ,(i/ B> «*) 2U^l(^a, ^b) F m (t)F; n ,(t) 

+ v/(m + l)(m' + l) 2iC 2 n 1 +1 , m , +1 K, vb) ^Ki^+iK, ^b) G m (t)G^(t), 

b*K,K,^;*) = uC("^) a t2-( l '^) f "(*ftW 

+ v /(m + l)(m' + l) 2 iC^ + i, m+ iK, ^b) SD^Uz+iK, ^b) G ro (t)G£,(t). 

In the next section we will discuss, in principle, three specific applications for the normally ordered moments with 
emphasis on: (i) the squeezing effect of the electromagnetic fields involved in this model; (ii) the nonlocal correlation 
measure which is based upon the total variance of a pair of EPR type operators; and finally, (iii) the Shchukin-Vogel 
criterion. 

III. APPLICATIONS 

Recently, Shchukin and Vogel [HI have further developed the concept of the complete characterization of single- 
mode nonclassicality (basically, this concept is based on the negativity of the Glauber-Sudarshan quasiprobability 
distribution function and requires an infinite hierarchy of conditions formulated either in terms of characteristic func- 
tions [ItJ or in terms of observable moments [Hj]) with the aim of characterising the entanglement of bipartite quantum 
states for physical systems described by continuous variables. In fact, they have provided necessary and sufficient 
conditions for the partial transposition of bipartite harmonic quantum states to be nonnegative, such conditions being 
formulated as an infinite series of inequalities for the moments of the states under investigation. Following the authors, 
'the violation of any inequality of this series is a sufficient condition for entanglement.' In this section, we will show 
that the moments derived for the driven JCM not only are considered as a theoretical realization of these conditions 
but also can be applied, for example, in the study of squeezing effect and in the numerical investigation of some recent 
proposals of inseparability criteria. 

A. Squeezing effect 

Let us introduce two dimensionless quadrature operators Q = (c 1 * + c)/v2 and P = i(c^ — c)/V2 for the 
electromagnetic field described by the q-number c, satisfying the commutation relation [Q,P] = il. In general, the 
squeezing effect occurs when either the variance Vq(£) = (Q 2 (i)) — (Q(t)) 2 or Vp(i) = (P 2 (i)) — (P(£)) 2 assumes 
values less than 1/2 - namely, Vcj(t) < 1/2 or Vp(i) < 1/2, but not both simultaneously (for more details on the 
squeezing effect and its applications in physics, see the review article in [43[). 

Figure 2 shows the plots of Vq (t) (solid line) and Vp' h \t) (dashed line) versus K ett t S [0, 50] when the atom-field 
system is resonant (a),(c) S — and nonresonant (b),(d) 5 = 6K cff , for e a = 3/\/l0, eb = 1/y/lO, \v a \ — 1, and 
|z/b| = 2 fixed. In particular, pictures (a) and (b) represent the variances associated with the quadrature operators 
Q a (i) and P a (*) (cavity field), while (c) and (d) correspond to the operators Q b (i) and Pb(t) (driving field). A first 
analysis of these pictures leads us to observe that the squee zing effect can be linked directly with the amplitude of 
the cavity (driving) field |^ a (b)| and the detuning frequency 5 |23l [il] . Indeed, numerical investigations for ^(b)! ^> 1 

and S > show that V^' b) (i) and Vp* ,b) (t) have the same asymptotic value 1/2, which minimizes the Heisenberg 

uncertainty relation VQ' h \t)Vp' h \t) > 1/4. Furthermore, an interesting aspect of this important effect refers to 

picture (d) , where the variance Vq (t) associated with the single-mode external quantum field exhibits the squeezing 
effect for S (from the experimental point of view, this theoretical prediction can be confirmed through balanced 
homodyne detections or similar techniques - for instance, see the experimental apparatus used by Josse and co-workers 
[Til in the demonstration of both quadrature and polarization entanglement generated via the interaction between a 
coherent linearly polarized field and cold atoms in a high finesse optical cavity). Beyond these important points on 
the squeezing effect, it is worth mentioning that the similarities between the different patterns of curves observed in 
(a),(c) for S — and (b),(d) when <5 = 6K off , can be explained by means of the entanglement dynamics verified in 
K ett t > among the constituent parts of the tripartite system (in this case, the continuous variable entanglement for 
the cavity and driving fields is obtained by tracing over the atomic variables). Thus, for the specific initial states of the 
cavity and driving fields adopted in this work, we can conclude that nonclassical effects and entanglement present an 
important link in the driven JCM, and this fact leads us to investigate some recent proposals of inseparability criteria 
for continuous bipartite quantum states. Next, we will focus upon the following question: "Can EPR uncertainty be 
considered the first quantitative characterization of the entanglement properties of the system under consideration?" 
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FIG. 2: Time evolution of variances associated with the quadrature operators Q a (b)(t) and P a (b)(t) for the cavity (driving) 
field considering the atom initially prepared in the excited state, and both the cavity and driving fields prepared in the coherent 
states. Pictures (a) and (b) represent the variances Vq^(t) (solid line) and Vp (i) (dashed line) for the following detuning 
frequencies: (a) S — (resonant) and (b) 5 = 6k c!! (nonresonant), with e a = 3/vTO, eb = l/\/T0, |i/ a | = 1, and \uh\ = 2 fixed. 
In addition, pictures (c) and (d) correspond to Vq (i) (solid line) and Vp (t) (dashed line) for the same set of parameters 
adopted in (a) and (b), respectively. The dot-dashed line present in all pictures refers to the variances of the coherent states 
(also described in literature as minimum-uncertainty states). Note that the variances show squeezing effect, this effect being 
connected with the amplitude |f a (b)| of the cavity (driving) field and the detuning frequency 5. 



B. EPR uncertainty 

For any bipartite physical system defined in a Hilbert space "K = 9£ a <g> 'Kb and described by an arbitrary 
normalized density operator p(t), the total variance of a pair of EPR-like operators - such as (Q a — Qt>)/\/2 and 
(P a + Pb)/V / 2, with Q a (b) and P a (b) satisfying the commutation relation [Q Q ,P^] = i5 a ^l for a,0 = a, b - allows 
us to establish the mathematical relation [y] 

I a ,b(f) = \ [v^\t)+V^\t)+V^\t)+vP{t)-2Cov^ h) (t) + 2Cov^' h \t)] >0 (12) 

where the covariance Cov^ b ' ) (t) for the operators X a g JC a and Xb € 1Kb is given by 

Cov^ b) (t) = ±Tr [p(t){Xa, X b }] - Tr [p(t)X a ] Tr [p(t)X b ] , (13) 

and Tr [p(i){X a , Xb}] represents the anticommutation relation mean value. Note that I a ,b(£) reaches the value zero 
when p(t) behaves as an idelized EPR type density operator. Furthermore, for states with Gaussian statistics, 
Ia,b(i) < 1 not only implies in the existence of nonlocal correlations but also establishes a sufficient condition for 
entanglement (in particular, this inseparability criterion has already been used in several experiments to demonstrate 
continuous variable entanglement - for instance, see [Llj and references therein). The additional condition I a ,b(i) > 1 
characterizes separable states into this context. Hence, the quantity I a ,b(i) measures the degree of nonlocal correlations 
and the inequality (fT2")) defines the EPR uncertainty for p{t). In Appendix A we determine explicitly the generalized 
moments (a^(£)a«(i)bt r (£)b s (t)), which are necessary in the numerical evaluation of the covariance function. 

Figure 3 shows the plots of I a ,b(i) versus K etf t for different values of detuning frequency: (a) 6 = 0, (b) S = 6k c!! , (c) 
S = 20K oft , and (d) S — 60ft off , where we have adopted the same set of parameters established in the previous figure. 
Note that EPR uncertainty is completely satisfied in all pictures, but unfortunately it does not reveal entanglement 
in the driven JCM since l&,b(t) > 1- According to Giedke and co-workers [l2|, the condition I a ,b(£) < 1 is met 'only if 
at least one of the uncertainties of (Q a — Qb)/v / 2 or (P a + Pb)/v / 2 lies below 1 (the standard quantum limit). This 
implies that the corresponding states must possess a certain squeezing.' In this sense, despite Vq^i) and Vq^t) have 
shown signatures of a nonclassical behaviour in pictures 2(b) and 2(d) for <5 = 6K otf , the inseparability criterion does 
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FIG. 3: Time evolution of I a ,b(£) versus K c!! t for e a = 3/%/TO, eb = l/\/l0, |i/ a | = 1, and \vt,\ =2 fixed, with different values 
of detuning frequency: (a) 8 = (resonant regime), (b) S — 6n ctt , (c) S = 20K cff , and (d) S = 60re e ff. Although the variances 
associated with the operators Q a (£) and Qb(i) have shown a nonclassical behaviour (squeezing effect) for the nonresonant 
regime, the EPR uncertainty can not be considered the first quantitative characterization of the entanglement properties in 
the driven JCM since I a ,b(£) > 1 for all t > (this case characterizes separable states). Beyond this important fact, it is 
worth mentioning that different initial conditions for the cavity and driving fields can be considered in the study of squeezing 
and entanglement. However, additional inseparability criteria must be employed in the quantitative characterization of the 
entanglement properties such as, for example, the Shchukin-Vogel criterion [la ]. 



not demonstrate the link between squeezing and entanglement in the theoretical model under study. Similar result 
was obtained by Shchukin and Vogel [16J for the example of an entangled quantum composed of two coherent states, 
whose expression reads as 



|0_)= 2 1-e 



-2(M 2 + |/3| 2 ) 



-1/2 



(\a,(3) 



a,-P)) 



In particular, the authors showed that the inseparability criteria established in [1,0] fail to demonstrate the entangle- 
ment of \ip-)- Consequently, our results lead us to investigate new inseparability criteria based in principle on certain 
series of inequalities involving higher-order moments which can be applied to a variety of quantum states. Thus, 
the next step will consist in the study of simple subdeterminants derived from the Shchukin-Vogel criterion [Til ] for 
entanglement. 



C. Shchukin-Vogel criterion 



Recently, Shchukin and Vogel [161 ] derived a hierarchy of necessary and sufficient conditions for the negativity of 
the partial transposition of bipartite quantum states, which are formulated in terms of observable moments associated 
with a variety of quantum states. This hierarchy is basically characterized by a set of inseparability criteria (or 
sufficient conditions for ent ang lement) that generalizes some previous purposes established in the literature (e.g., see 
Refs. 0,0!! EM, EE [l9|). For practical reasons, let us formulate the Shchukin-Vogel criterion as follows: 
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FIG. 4: Time evolution of S(t) versus n elt t for different values of detuning frequency: (a) 5 = (resonant regime) and (b) 
S = 6re e tf (nonresonant regime), where we have adopted the same set of parameters established in Figure 2. Note that S(t) 
assumes negative values in both pictures and this fact corroborates the Shchukin-Vogel criterion for entanglement. 



Shchukin-Vogel criterion. The partial transposition of any bipartite quantum 
state is nonnegative if all the determinants derived from the Nth order determinant 



1 
(at) 
(a) 
(b) 
(bt) 



r N = 



(a) 


(at) 


(bt) 


(b) 


ata) 


(af> 


(atbt) 


(atb) 


(a 2 ) 


(aat) 


(abt) 


(ab) 


(ab) 


(atb) 


(btb) 


(b 2 ) 


abt) 


(atbt) 


(b^) 


(bbt) 



are nonnegative (in other words, VN we obtain T N > 0). However, if there exists 
a negative determinant (namely, 3N such that T N < 0) the negativity of the partial 
transposition has been demonstrated - this fact provides a sufficient condition for any 
bipartite quantum state under investigation to be entangled. 

To illustrate this criterion within the driven JCM we consider, in principle, the time-dependent subdeterminant 



S(t) 



1 

(b(*)) 
it(i)b(t)) 



(bt(*)) 
(bt(i)b(t)) 
it(t)bt(i)b(t)) 



(a(t)bt(t)) 
(a(*)bt(t)b(i)) 
it(t)a(t)bt(t)b(*)) 



(14) 



where the mean values (at p (t)a 9 (t)bt r (f)b s (t)) can be evaluated through the results obtained in Appendix A. Figure 
4 shows the plots of S(t) versus K et! t for e a = 3/V10, £b = 1 / a/10, |^a| = 1, and |z/b| = 2 fixed, with different values 
of detuning frequency: (a) (5 = (resonant regime) and (b) S = 6ft off (nonresonant regime). From these pictures 
we can observe that S(0) = (this condition refers to separability of the initial density operator for the tripartite 
system) and when n o{1 t > 0, the time-dependent subdeterminant l|14p assumes negative values in both situations (the 
entanglement condition is promptly verified). In particular, the inseparability of the reduced density operator turns 
to be more evident in the nonresonant regime - see picture 4(b) - which coincides with the squeezing effect observed 

in pictures 2(b) and 2(d) for the variances Vq (t) and Vq (t), respectively. Indeed, the link between squeezing effect 
and entanglement turns more evident through the detuning frequency S = ujq — to, this frequency being associated 
with the atomic transition (coq) and cavity (driving) field (u) frequencies in the driven JCM. Thus, by means of 
the numerical evidence obtained in this section, we can conclude that S(t) corroborates the Shchukin-Vogel criterion 
and consequently, it can be considered the first quantitative characterization of the entanglement properties for the 
theoretical model under investigation. 



IV. CONCLUSIONS 



In this paper we have applied the decomposition formula for su(2) Lie algebra on the driven Jaynes-Cummings 
model in order to calculate the exact analytical expressions for the normally ordered moments when the atom is initially 
prepared in the excited state. In particular, adopting the diagonal representation of coherent states, we have shown 
that the Wigner characteristic function can be written in the integral form, with their integrands having a commom 
term which describes the product of the Glauber-Sudarshan quasiprobability distribution functions for each field, and 
a kernel responsible for the entanglement. Next, we have evaluated the normally ordered moments for the cavity and 
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driving fields separately, by means of their respective Wigner characteristic functions. It is worth emphasizing that 
the mathematical procedure adopted here does not present any restrictions on the states of the cavity and external 
electromagnetic fields. To illustrate our exact results we have fixed, for convenience in the numerical calculations, 
both the cavity and driving fields in the coherent states. This procedure has allowed us to investigate not only the 
squeezing effect associated with the quadrature operators Q a r b ) (t) and P a (b) (t) of the cavity (driving) electromagnetic 
field, but also the EPR uncertainty and the Shchukin-Vogel criterion. Such numerical investigation has produced some 
important results, within which some deserve to be mentioned: (i) we have shown that squeezing effect is directly 
related to detuning frequency in the driven JCM; (ii) the EPR uncertainty has failed in the description of entanglement 
properties of the system under consideration; and finally, (iii) we have demonstrated that the Shchukin-Vogel criterion 
correctly establishes a set of necessary and sufficient conditions for the negativity of partial transposition of bipartite 
quantum states. 
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APPENDIX A: EXACT ANALYTICAL EXPRESSIONS FOR THE GENERALIZED MOMENTS 



The main aim of this Appendix is to obtain an exact analytical expression for the generalized moments 

(a tp (i)a«(t)b tr (i)b s (i)) ee Tr [p(t) a tp a 9 b tr b s ] (Al) 

with {p, q, r, s} £ N, by means of the Wigner characteristic function 

AwKa.C&.&t) = Tr[p(i)D a (£ a )D b (£ b )] (A2) 

which involves the expectation value of the product of displacement operators for each field in the driven JCM. To 
this end, let us establish the link between both expressions through the mathematical relation [4l| 



(at"(t)a«(*)bt'-(i)b a (t)) = (-1)<H 



QP+q+r+s 



■ exp 



^(ICaP + ICbl 2 ) 



Aw(£a> £ a > £b) £b! 



(A3) 



Thus, our first task will consist in obtaining the Wigner characteristic function (|A2[) . while the second task will focus 
upon the evaluation of the generalized moments (|Alj) with the help of equation (IA3|) . 

Adopting the same mathematical recipe used in Sec. II to derive Eq. ([3]), the generalized Wigner characteristic 
function is expressed into this context as 



A w (£a, = // d ^ p a(a a )Pb(ab)3Cg a(b) ,Q (a ai a h ;t) , 



(A4) 

where P a (a a ) and P b (a b ) are the Glauber-Sudarshan quasiprobability distribution functions associated with the cavity 
(a) and driving (b) fields, and 



X 



«a ( b),Q b) K,0!b;*) = nD| a 1 (b)iC(b) (a a ,a b ;i)+ 2iD£ (b)lC; K, a b ;t) 



(A5) 



represents a kernel responsible for the entanglement between the displacement operators D a (£ a ) and D b (£ b ) with 

ijD" „ K,a b ;*) = (a a ,a b |Ul(t)D a (e a )D b (a)U ij (t)|a a ,a b ) . (A6) 

S£ H b )' s a(b) J 

After lengthy calculations, the analytical expressions for the mean values (|A6|) assume exact forms similar to that 
obtained for the cavity field but differ in the dependence on the variables £ a ( b ) and £ a ( b )- As a consequence of this 

fact, the complex function Y^"^ m \a a , a b ) must be adequately substituted by 

^"^Tf* Ob) = exp { -- (|£ a | 2 + |£ b | 2 ) + 2ilm [(e b a a - e a ab)* (fb^a - e a Cb 



x [(e a a a + e b a b )* (e a £ a + e h £ h )] m ' m L^'^ (|e a ^ a + e b ^ b | 2 ) 
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for each situation described by the indices i,j= 1,2. In particular, if one considers the initial states of the cavity 
and driving fields in the coherent states, we obtain A w (£ a , £ a , £b, £ b ; t) = ~^£a(b) >£J (b) ( v *i an< ^ consequently this 
characteristic function assumes the simple form 



(A7) 



Following, we will use (|A4|) to derive an exact analytical expression for the generalized moments (jA3|) . 
Let us initially introduce the generating function for < K^ n ' m X (a a ,a b ) through the definition 

Sa(b) >£ a (b) 



1 -^m+m 

Sa(b)'«a(b) TO! 9w m 3t>" 1 ^"( b ) ^a(b) x 



li.D— 



where 



u k v k 



k,k'=0 



Furthermore, we also consider the auxiliar function 



p.p+q+r+s 

«(P,Q,r,s), \ _ (_y\q+s " 



(l^aP + ICb 



nj-im,m ) i \ 



(A8) 



(A9) 



(A10) 



which has a fundamental role in the present context since this function should reduce the complexity of our calculations. 
Indeed, substituting (|A8| into Eq. (|A10p we obtain promptly the following useful relation: 



Ira ni* K, <*b) 



1 3 



m+m 



i! bu m dv' 1 



^afe«^ s) K,«b) 



(All) 







with du,v' V ' S) {ota., a h ) given by 



a^)(a a ,a b ) = (-l) 



g+s _ 



gp+9 



+r+s 



^a| 2 + |£b| 2 ) 



n (u,v) I \ 

«a(b) 'Sa(b) 



«a,«*,«b,«*=0 



For instance, the generating function (|A9[) reads as 

s Sbte; (b ) (aa ' ab) = exp {~l ^ a|2 + l?b|2 ) + 2iIm [ (£btta _ £aC * b) * (£b?a " £a&) ] } 

(e a £a + e b £ b )* 



(A12) 



(e a a a + e b a b )* 



x exp 

while (|A12[) has the explicit expression 

du,v' r ' s) ( a ^ a h) = [e a (e a a a + £b«b)* v + e b (e b a a - e a a b )*] P 
x [eb (e a a a + e b a b )* v - e a (e b a a - e a a b )* 



u + (e a a a + e b a b )* (e a £ a + e b £ b ) v + uv 



ea" 



e a a a + eb"b) 
e b u 



(e a a a + eb«b) 



+ £b (£b" a - e a "b 
e a (eba a - e a ab) 



exp(itv) 



Consequently, the partial derivatives of this expression with respect to u and v at the point u, v — allow us to 
determine the auxiliar function (|A11[) . 

Finally, it is worth mentioning that the generalized moments for any kind of cavity and driving fields can be obtained 
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by means of a basic sequence of functions evaluated in this appendix - e.g., the schematic diagram showed below 



l J t»(b)i« a(b ) ^ a ' 

^L,b),Q b) K,ab;i) 

A w (£a, ^a; Cbj Cbi t) 



nj-(m,m') , \ 



St"'" e. ((*a,Qb) 

Sa(b) iS a ( b ) 



T 
I 

Y 



du,v (Q!a,Q!b) 



T 
I 

Y 



n (p,q,r,s) / \ 
-A 



ltP(£) a 9(t)bt r (t)b s (*)) 



represents two alternative ways of finding out the generalized moments (|A1[) ■ From the theoretical point of view, (|A1[) 
plays an essential role not only in the investigation process on both the squeezing effect and the EPR uncertainty, 
but also in the analysis of some inseparability criteria for the driven JCM. 
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